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Abstract 

We discuss physical and mathematical aspects of the over-damped motion of a Brownian particle 
in fluctuating potentials. It is shown that such a system can be described quantitatively by fluctuating 
rates if the potential fluctuations are slow compared to relaxation within the minima of the potential, 
and if the position of the minima does not fluctuate. Effective rates can be calculated; they describe 
the long-time dynamics of the system. Furthermore, we show the existence of a stationary solution 
of the Fokker-Planck equation that describes the motion within the fluctuating potential under some 
general conditions. We also show that a stationary solution of the rate equations with fluctuating rates 
exists. 

1 Introduction 

Thermally activated relaxation processes are an important mechanism for the dynamics of many physical, 
chemical, and biological systems. The diffusion of a Brownian particle in a potential with several minima 
and barriers serves as a paradigm for such relaxation processes. The minima of the potential represent 
stable or metastable states of the system, and the motion of the particle models the transition of the system 
from one state to the other and back. The most simple example in this class of models is the problem 
of diffusion over a single potential barrier, pioneered by Kramers [|l|]. The dynamics of the diffusion 
over a barrier is dominated by a characteristic time scale, which is given by the mean first passage time 
for the escape out of the minimum of the potential. In a model with several minima and barriers, the 
corresponding time scales are given by the inverse rates for the transition from one state to another. 

Often one uses kinetic rate equations to describe the dynamical properties of relaxation in chemical or 
biological systems. This description is clearly much simpler than the description by a Brownian particle 
in a potential. To calculate rate coefficients, it is often sufficient to know some of the properties of the 
potential energy surface, like barrier heigths. In principle, if one knows the potential energy surface, both 
descriptions are equivalent. One can calculate static or dynamic properties of the system either using the 
Brownian particle in the potential or using the rate equations. 

In many situations, the potential fluctuates due to some external fluctuations, chemical reactions, or 
oscillations. These fluctuations usually have a finite correlation time, which corresponds to a non-thermal 
noise. The most simple model with a fluctuating potential consists of a single, fluctuating barrier, where 
the height of the barrier fluctuates between two different values. Doering and Gadoua investigated 
such a simple model. They found a local minimum in the mean first passage time as a function of the 
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barrier fluctuation rate. This effect has been called resonant activation and has been studied extensively 

In principle it should be possible to use a description with rate equations in the case with fluctuating 
barriers too. This picture has already been suggested by Bier and Astumian [Q]. They showed that the 
long-time dynamics of a simple model with a dichotomously fluctuating linear ramp can be described 
by kinetic rate equations if the potential fluctuations are not too fast. A similar observation was made in 

where it was shown that noise induced stability of a metastable state occurs in systems with a single 
fluctuating barrier. 

In the present paper we consider potentials in d dimensions with several minima and barriers. Fur- 
thermore, we allow general potential fluctuations. For simplicity we assume that the potential fluctuations 
can be parametrized by a single Markovian noise process z{t). Often it is argued that the Markovian noise 
z{t) represents the cumulative effects of many weakly coupled environmental degrees of freedom. In that 
case the central limit theorem can be applied and z{t) becomes an Ornstein-Uhlenbeck process. On the 
other hand, there are many realistic situations where the fluctuation of the potential is triggered by a single 
(or few) environmental degree of freedom. The most simple case of a non-Gaussian fluctuating potential 
is a dichotomously fluctuating potential as discussed e.g. in But in many applications, the potential 
fluctuates not only between two different values. Since all these different cases are of interest, we will 
investigate general noise processes. Furthermore, it is clear that there are always many possibilities to 
parametrize a fluctuating potential by a Markovian noise process. This is another reason why it is useful 
to consider general noise processes. 

A main purpose of the present paper is to show under which conditions it is possible to pass from 
a model with a fluctuating potential to a set of rate equations with fluctuating rates. We will discuss 
general models in arbitrary dimensions and with potentials having several minima. The general picture 
we support by our calculations is that rate equations can be used if the potential fluctuations are slow 
compared to the relaxation within in minima of the potential and if the positions of the minima of the 
potential do not fluctuate. The long-time behaviour of the model with a fluctuating potential is then 
similar to that obtained using rate equations. The fluctuation of the potential is mainly a fluctuation of the 
barrier heights. The potential fluctuations may be fast compared to relaxation across the barriers of the 
potential. 

After these remarks it is clear but important to mention that the applicability of rate equations depends 
also on the quantities one in interested in. As long as one is interested in stationary or quasi-stationary 
properties, rate equations yield good results. If one is interested in properties on time scales of the order 
of the relaxation within the minima or smaller, rate equations cannot be applied. 

Furthermore we analyse systems with fluctuating rates in detail and derive effective rates (which do 
not fluctuate). The effective rates yield the relevant time scales for the long-time behaviour of the system. 
They can be directly compaied to the mean escape rates in the fluctuating potentials. Furthermore they 
yield the stationary distribution, which can also be compared with results obtained directly by solving 
the Fokker-Planck equation with a fluctuating potential. The comparison is easily done for simple, one- 
dimensional systems with piecewise linear potentials. For such systems it is relatively easy to solve the 
Fokker-Planck equation directly. 

The paper is organized as follows: In the next section we discuss how a kinetic rate equation with 
fluctuating rates can be obtained starting from a Fokker-Planck equation with a fluctuating potential and 
in which situations one can expect it to be valid. We also show how one can derive effective rates from 
the kinetic rate equation with fluctuating rates. The effective rates yield the chaiacteristic time scales for 
the escape out of the minima of the fluctuating potential. The general formula for the effective rates can 
be evaluated for special noise processes. This is done in section 3. The results are used to compare the 
description of a system by effective rates with the original Fokker-Planck equation. For our discussion 
the existence of a stationary solution of the Fokker-Planck equation with a fluctuating potential and the 
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existence of a stationary solution of a kinetic rate equation with fluctuating rates is essential. Section 3 
also contains the exact solution of the dichotomous two-state system and its derivation. Some of the exact 
results have already been mentioned without derivation in [^]. In section 4 we show that the stationary 
solutions of the rate equation with fluctuating rates and of the Fokker-Planck equation with a fluctuating 
potential exists under some general assumptions on the potential. This result is a generalization of a 
recent, similar result by Pechukas and Ankerhold [|l^]. Finally, in section 5, some possible applications 
are mentioned and conclusions of our results are drawn. As a specific example we discuss how our results 
can be used to interprete experimental findings for membrane proteins. 

2 Derivation of effective rate equations 
2.1 Noise processes and the Fokker-Planck equation 

The over-damped motion of a particle in a fluctuating potential is usually described by a Langevin equa- 
tion 

f = /(x,o + V2r|(0- (1) 

at 

We have chosen the units such that the friction constant is unity, ^{t) is a thermal (Gaussian white) noise, 
it satisfies = 0, = 5„/,5(f - 1'). f{x,t) = -VV(x,f) is the force of the fluctuating po- 

tential. We investigate situations where the potential tends to infinity for \x\ ^ oo or where the potential is 
defined in a finite domain ^l with reflecting boundary conditions on the boundary dQ.. Let us assume that 
the potential fluctuations can be parameterized by a single stochastic variable z{t), V{x,t) = V{x,z{t)). It 
is clear that such a representation is not unique. Instead of z(f) one can use any monotonic function of z{t) 
to parametrize V{x,t). Often one uses a simple linear ansatz V{x,z{t)) = Vq{x) + zit)AV (x), with a di- 
chotomous or an Omstein-Uhlenbeck process z{t). As already mentioned in the introduction, we want to 
investigate models for which the potential fluctuations represent the effects of few environmental degrees 
of freedom. Therefore we dot not restrict ourselves to a linear ansatz or to special noise processes. The 
only assumption we make is that the noise process z{t) can be described using a Fokker-Planck equation 
for the probability distribution p{z,t), 

'-^=MMzj). (2) 
The right eigenfunctions of the generator M-. are denoted by (|)„ (z) , the eigenvalues by — ?i„ 

MMz) = -K<^n{z), (3) 

where Xq = 0, Xi > 0, X„ > (|)o(z) is the stationary distribution of the noise process z(f). We let 

<^niz) = g„(z)(t)o(z)- Then one has 

j dzgn{z)gm{z)^i){z) = K,m- (4) 

gn{z) are orthogonal functions with respect to the weight function (|)o(z)- For n = m this equation fixes 
the normalization of (|)„(z). The eigenvalue X\ = x^^ determines the correlation time of the noise process 
z{t). We already mentioned that the pai^ametrization of the fluctuating potential by a stochastic variable 
z{t) is not unique. One can always choose a different parametrization. The final results should of course 
be independent of the parametrization. If one changes the parametrization, the spectrum of remains 
the same, but the right eigenfunctions are changed. 
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The motion of the over-damped particle in the fluctuating potential can be described using the joint 
probabiUty distribution p{x,z,t) for and z{t). It obeys the Fokker-Planck equation 



' ' = V • {{yV{x,z)) + TV)p{x,z,t) +M,p{x,z,t). (5) 

V{x,z) must have the same properties as a function of x as V{x,t). It must tend to infinity sufficiently 

fast when |x| becomes large; or alternatively, we may define the Fokker-Planck equation on a finite, open 
domain Q. with reflecting boundary conditions on the boundary dO.. Furthermore we assume that V{x,z) 
is finite for any z out of the support of <^o{z). V{x,z) does not have any infinitely high potential wells 
inside Q.. 



2.2 Rate equations 

2.2.1 Rate equations for fixed potentials 

Let us first briefly review the case of a potential that does not fluctuate. The Fokker-Planck operator has 
the form 

L = V-((VV(x)) + rV). (6) 

Let us assume that the potential V{x) has A'^ distinct minima at the points x,, / = 0, . . . ,A'^ — 1 and that 
V{x) — > oo for 1^1 — > oo sufficiently fast. L has one eigenvalue 0, the corresponding right eigenfunction 
is the stationary distribution <x exp{—V{x)/T). — L is non-negative, i.e. all the other eigenvalues of —L 
are positive. Using a simple variational argument, one can show that —L has at least (A'^ — 1 ) eigenvalues 
which behave like exp(— c/r) for small T. Let us divide the entire space on which the system is defined 
into subspaces xi G £lj if and only if / = j. Q.i n ily = if / / j. U/^i = \ {{Ji^^i}- denotes 
the boundary of H,. The boundaries are chosen such that • VV{x) = for x G 3£2,-, where is 
the normal vector to the boundary at x. Such a division of is always possible. Depending on the 
form of V (x) it may not be unique. Let us now introduce subsets Cli C il, which do not contain a small 
region close to the boundary of but which do contain the positions Xi of the minima together with 
a sufficiently large region around them. Let 9, (x) be a smooth function that is equal to unity if x G H,- 
and that vanishes if x ^ Q.,. In il, \ Cl, the function e,(x) drops smoothly from 1 to 0. We assume that 
sufficiently many derivatives of 6, exist. 

In a next step, we introduce L = exp(V(x) / (2r)Lexp(— V(x) / (27)). L is a hermitian operator. It has 
the expUcit form 

L = rv2 + l(vV)-^(vy)2. (7) 

Clearly, L has the same eigenvalues as L. exp(— V(x) / (27)) is the (unnormalized) eigenfunction of L to 
eigenvalue 0. Let g,(x) = C,exp(— V(x)/(2r))6,(x), where Q is chosen such that / d'^x\gi{x)\^ = 1. One 
has / d'^xgi{x)gj{x) = dij. Let Lij = J d'^xgi{x)Lgj{x). These quantities vanish for i ^ j. For i = j one 
obtains 

Ui = I IQp I d^xv . (exp(-y (x)/r) ve,-(x)). (8) 

Since V6,- vanishes forx G Cli, Lu contains a factor exp(— At^/T) where 

AVi= min V{x)-V{xi). (9) 

x€Cii\d.i 

Let Ay = min, AV/. Then this shows that —L has eigenvalues that decay at least as fast as exp(— AV/r) 
to zero when T tends to zero. 
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Let us now consider L restricted to but with reflecting boundary conditions on 3n,. With these 
boundary conditions —L has one eigenvalue and typically other eigenvalues which are much larger than 
exp(— AV/r). These eigenvalues describe the time-scale for a decay within the potential well around 
the minimum at Xi. On the original domain Q. one can therefore expect that — L has one eigenvalue 
and A^— 1 eigenvalues which are of the order exp(— c/T) and that all other eigenvalues of —L are much 
larger if T is sufficiently small. Let us denote the smallest eigenvalues of —L by A.,, ? = 0, . . . , — 1, 
Xo = and the corresponding right eigenf unctions by gi{x). Then, for time scales large compared to the 
intra- well relaxation times of the potential we can make the ansatz 

N-\ 

P(^>0 = L 9i{t)gi{x) (10) 

(=0 

for the solution of the Fokker-Planck equation ^ = Lp. This ansatz should describe the long-time be- 
haviour of the system quite well. It is valid on time scales being large compared to the intra-well re- 
laxation times in the different potential wells of the potential. = ^^^Q^,d'^xp{x,t) is the probability 
to find the particle in the region il,-, i.e. close to the minimum Xi. Let Uij = SxeQ.i^'^^Sii^)- has 
n,(f) = Y^j Pj{t)nij. For the derivative with respect to time one obtains 

^ = f d'xLp = Y,pj[ d'xLgj 

at JxeQ-i j JxeCii 

= ~ L ^iP = Y^nknk, (11) 

j k 

where the matrix R = (?',j),j=o,...,A?-i is given by R = —N^^AN. N is the matrix («(7);,;=o,...,Af-i and A is 
a diagonal matrix with the entries Xj, j = 0,...,N—\. The differential equation ^ = Y,k ^ik^k is the rate 
equation for this problem. By construction, the eigenvalues of R are —X,. 

Unfortunately, the definition of R is not useful for an explicit calculation. If one is able to solve the 
eigenvalue-problem of L, there is no need to describe the problem by a rate equation. The rate equation 
is useful in situations, where it is not possible to obtain analytic results from the Fokker-Planck equation. 
Fortunately, approximative values for the rates r/y often yield very good results. In the one-dimensional 
case the usual Kramers' rate is sufficiently good. 

By construction, R is equivalent to the Fokker-Planck operator L projected onto the space of its eigen- 
functions with largest eigenvalues. This means that dynamical properties on time scales corresponding 
to eigenvalues or larger can be well described by R. Properties on time scales IX^rl"^ or smaller cannot be 
calculated using R. Xn sets the time scale for relaxation within the potential wells. 



2.2.2 Rate equations for fluctuating potentials 

Under which conditions is it possible to describe a problem with a fluctuating potential by rate equations? 
It is clear that one must be able to define n,- for a fluctuating potential. This is possible if the stationary 
distribution is strongly peaked at the minima x,. But in general one can have a situation, where the 
positions x, of the minima fluctuate as well. If this is the case, the probability distribution will not 
necessarily be peaked. Let us discuss that using a simple example: A one dimensional dichotomously 
fluctuating potential. 

For this example, z takes two values, z = ib 1 with equal probabiUty. has the eigenvalues and 
— T~\ In the eigenbasis of the potential V{x,z) has the form 

( {V){x) AV{x) \ 

{AVix) {V)ix)J^ ^'^^ 
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where {V) {x) = {V{x,+) +Vix-))/2, AV{x) = -y(x,-))/2. Let /(x) = A/(x) = 

~^^x^- calculate the stationary probability poix) = J dzp{x,z). The stationary joint probability 

distribution can be written as p(x,z) = Po{x)^o{z) + Pi {x)^i (z). Inserting this into the stationary Fokker- 
Planck equation yields 

^{f + zAf- T-^){po{x)Uz)+Pi{x)h{z)) +T-Vi = 0. (13) 

This yields two equations for po{x) and p\(x) corresponding to the two coefficients of (|)()(z) and (|)i(z), 
which both have to vanish. Ehminating pi{x) yields a single equation for po{x), which can be written as 



dx \ dx 



(A/)-' (/-r^jpoW + x^A/poW=0. (14) 



This equation can be derived directly from the Fokker-Planck equation. For T = it can be solved 
explicitely. One obtains 



^0^=^723^72^^? {-I J ^^723X72 1 • (15) 



The integrand // (/^ — A/^) diverges at the minima of V{x, it). If the position of the minima does not 
change, pq{x) is a sum of 6-function located at the minima with appropriate weights. With changing 
positions of the minima, the situation becomes different. In that case the particle moves between the 
minima. The most simple example to illustrate that is V{x,±) = \y±{x — x±)'^, y± > 0. Letx+ > x_. 
Then the support of po{x) is the interval For po{x) one obtains 

= ' (^^5^ + uhr^) ('-')"'^''*'(-.-)"'^-' (16) 

with some normalization constant C. The distribution function po{x) is not peaked at x±. It is clear that 
for small T > one obtains a similar result. This simple example shows that in a situation where the 
position of the minima depends on z the stationary distribution of the system is not necessarily peaked at 
the minima of the potential. 

Let us now assume that the potential V (x, z) has minima at Xi and that the positions Xi of the minima 
do not depend on z. For a fixed potential, rate equations can be used for time scales larger than the intra- 
well relaxation times. It is therefore clear that for potential fluctuations being faster than the intra-well 
relaxation times one cannot use rate equations. If the potential fluctuations are slower than the intra-well 
relaxation times, and if the position of the minima does not fluctuate, one is in an adiabatic situation. The 
particle is able to follow the fluctuations of the potential as long as it stays close to a minimum. This 
means that one can calculate the rates for fixed z- The rates become functions of z. The system is thus 
described by rate equations with fluctuating rates, 

^=LnMt))nj. (17) 

These rate equations are considerably simpler than the Fokker-Planck equation with a fluctuating po- 
tential. On the other hand, the validity of these equations is not entirely clear. We used only heuristic 
arguments to obtain the rate equations. In the following section we will compare results obtained from 
the rate equations with results from the Fokker-Planck equation. This allows us to show, in which range 
of parameters the rate equations can be used. 
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2.3 Elimination of 

If one is interested in stationary properties or in the long-time behaviour of the system, it is useful to 
derive effective rates for the transition from one state to another. The effective rates do not depend on 
z. They can be calculated by eliminating z{t). If the correlation time x of z{t) is small, the elimination 
process is a standard elimination of a fast variable. To lowest order in x the effective rates are simply 
given by the average of the rates r;;(z). To derive effective rates, let us transform the rate equation 

dn ^ 

— (18) 
at 

into a Fokker-Planck equation for the joint probability density p(n,z,t) 

dp{n,z,t) 



dt 



yn-{Rnp{n,z,t))+M,p{n,z,t). (19) 



Since = 1 for all t, p{n,z,t) contains a factor 8(^,?2,- — 1). To calculate average values for n, it is 
sufficient to consider the stationary case. Let 

p{n,z) = po{n)i^oiz)+pin,z) (20) 

where 

ldzp{n,z)=0, (21) 

J d^np{n,z)=0. (22) 
Let {R) = J dzR{z)<^o{z). The stationary Fokker-Planck equation yields 

Vh • ((/?) fipo{n)) = -Vn J dzMp{n,z) (23) 

and 

p{n,z) = -M;^VfiRn{poin)<^o{z)+p{n,z)). (24) 

Here, is the generalized inverse, it obeys M^^(^„{z) = — X^^(|)„(z) for n > 0, M^^(^o{z) = 0. (23) 
yields 

{R)npo{n) = - J dzRnp{n,z) + jin) (25) 
with Vnj{n) = 0. Let (n) = J d^nnpo{n). Since 

|j^n;(H) = ^^^=0 (26) 

one has 

{R){n) = fdz fd''ni-l)'Rn{M;'VnRny+'poin)i^oiz) 

= -f, f dz f d'^niRM^'y+'Rnpomoiz)- (27) 
1=0-' 



The first line in ( |27D is obtained by solving ( |24| ) iteratively and inserting the result in (|25|). Integrating by 
parts yields the second line. Since 7? and Mj do not depend on n, the integration over n yields the average 
(n). The sum over / gives finally 

J dz{l-RM-Y^R<^oiz){n)=0. (28) 
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The entries of the matrix 



Rcff = J dz{\-RM-^) ^/?(^o(z) 



(29) 



are the effective rates, the negative inverse of its diagonal matrix elements are the mean escape times out 



of the states /. The expansion in (g7|) can be viewed as a usual x-expansion. Expanding the right-hand 
side of ( [29| ) yields the x-expansion for the effective rates. But since (1 — RM^^) is positive definite, the 
right hand side of ( |2^ is well defined for any x. Whereas it is not possible to speak of an effective barrier 
for large x [|5|], it is possible to obtain effective rates. But, depending on the noise process and on R{z), 
it cannot be excluded that the smallest non- vanishing eigenvalue of /?eff tends to zero in the limit x — > oo. 
This may happen if the fluctuations of the barrier heights in the potential can be arbitrarily large, so that 
infinite barrier heights occur. In that case the solution for (H) has to be obtained from the solution for 
finite X and the limit x ^ °° has to be taken afterwards, since (28) does not have a unique solution for 
X = oo. As long as the eigenvalue of R{z) is non-degenerate for all z in the support of ^o{z), the matrix 



/?eff has a non-degenerate eigenvalue and (28) has a unique solution. This is shown below in section 4. 



2.4 Applicability of effective rate equations 

As we already discussed in subsection 2.1, rate equations describe only properties on long time scales. 
If one wants to investigate dynamical properties on characteristic time scales of the order of the typical 
relaxation within a potential well, one cannot use rate equations. If one wants to study a system with a 
fluctuating potential that fluctuates fast (i.e. one a time scale of the order of the typical relaxation within 
a potential well) one has to use the Fokker-Planck equation. This is true for equations with fluctuating 
rates as well as for the effective rate equations we derived in the last subsection. 

Since in many applications details of the potential are unknown, one often uses rate equations. Es- 
pecially in biological applications, e.g. for the description of membrane proteins, one usually uses rate 
equations. A typical example is the well known paper by Petracchi et al [|ll|], who studied the effect of 
time dependent electric fields on membrane proteins experimentally and used rate equations to interprete 
their experimental findings. They were able to explain most of their results using rate equations and 
a very simple ansatz for the rates. In their experimental studies they found an interesting effect called 
phase anticipation. The rate equations on the other hand do not show this effect. They argued that phase 
anticipation occurs due to some special biological effect. From our discussion about the validity of rate 
equations it is clear that the effect of phase anticipation, which occurs on short time scales, cannot be 
described by rate equations. We will come back to this point in our conclusions. 

In the following section we discuss some simple, one-dimensional examples. The most simple case, 
the dichotomous two-state model, is exactly solvable. For this case and for more general examples we 
compare the solution of the rate equations with the solution of the Fokker-Planck equation. For all cases, 
rate equations yield very good results if the potential fluctuations are slow compared to the relaxation 
within the potential wells. 



3 Solutions for special cases 
3.1 The dichotomous two-state model 

The most simple model with a fluctuating potential is a dichotomous two-state model. The potential 
V(x,z) has two minima at x\ and X2. z{t) is a dichotomous process. Such a model has two states, and 
corresponding rates r,y, /,7 = 1,2 for transitions between them . For this model we let ri := r2\ = —rw, 
ri '■= r\2 = —r22- The rate equations are 
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dni 

__ = -nm +r2?i2, (30) 
at 

__ = -r2?i2 + n«i- (31) 
at 

Using the normalization condition +?i2 = 1 one obtains a single differential equation 

^ = -{ri+r2)ni+r2. (32) 

In the following the index of «i will be dropped. The Fokker-Planck equation for p{n,z,t) is 

dp(n,z,t) d , s , N , 

^) = ^{{ri+r2)n-r2)p{n,z,t)+M,p{n,z,t). (33) 
at an 

In general it is possible to expand the stationary solution using the eigenfunctions of Mj 

p{n,z)='£PkinMz). (34) 
k 

For the matrix elements of r, (z) one has the representation 

riiz)<\>kiz)=Y.mUz)- (35) 
/ 

It is useful to define 

ru ■=nH + r2ki- (36) 
The Fokker-Planck equation for the stationary distribution p{n,z) has now the form 

a d 

ki 

Y.hPi<^i. (37) 



^{iri+r2)n- r2)p{n,z) = ^T,{rkin-r2ki)pk^l 
an an 



The two matrices Rt = {rju) have the same eigenvectors and only positive eigenvalues. This must be true 
for /? = /?i +7?2 as well. In the dichotomous case R and /?2 are 2 x 2-matrices. The eigenvalues are r± 
and r2±. ri{z) takes the two values r,± with the probabilities p±. Let n± = r2±/r±. I assume n^ > n^. 
The explicite form is of R is 

^^^roo m\( p^+p^r ./P^{r^-r^)\ 
V no rn J \ y/pTp^{r+-r^) p+r^+p^r+ J 



and similarly for /?2- Poin) vanishes for n <n- and for n> n+. The two equations ( p7| ) are 

{roon - r2oo)po + {rmn - r2w)Pi = 0, (39) 

{{roin - r2Qi)pQ + {run - r2n)pi) = 'k\p\. (40) 

an 

The integration constant in ( 09| ) vanishes since the average of ^ vanishes. It may be used to ehminate 
p\, one obtains 

^ApQ = -hBpQ (41) 
an 
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where 



no" -''210' 



(42) 



no" - '•210 



The solution of this equation is 



770 = -exp(->,i j -dn) 
where C is a normalization constant and f{n) is given by 



(44) 



f{n) := (roi?i-r2oi)(rio?i-r2io) - (roo«-r2oo)(ni«-'"2ii) 
= -det{Rn-R2) 

= -det{R){n-n+){n-n^). (45) 



The integral in (44) can be calculated. The explicit solution is 

Po{n) =Cin-n){n-n_f--\n+-n)''^-^ (46) 

where 

a_ = ^^'^f^-^-\ , (47) 
= (48) 



r+r_(n_|_ — «-) 
''210 r2+ - r2- 



(49) 



One can show that either n <n^ or h> n^, so that pQ{n) is non-negative for n € as it should be. 

For Xi —>■ 00 one obtains the expected result Po{n) = 6(« — «o) where no = r2oo/''oo- For small Xi, po{n) 
has algebraic singularities at the boundaries n±. For small but finite Xi , the singularities remain integrable. 
In the limit Xi —>■ 0, po (n) tends to a sum of two 5-functions (with appropriate weights according to the 
dichotomous process), located at n± = r2±/r±, 

PQ{n)=p+?>{n-nj^)+p^^{n-n^). (50) 

The average value hofn can be calculated explicitely. One obtains 

npQ{n)dn 

n{n -n){n- - n)''^-^dn 

f^^ {n — n){n — n^)'^-^'^ {n+ — n)"-+^^dn 

n^{n^ — «)B(a+,a_) + (2w_ - n){n+ — «_)B(a+,a_ + 1) + («+ — «-)^B(a+,a- +2) 
(«_ - n)B{a+ ,«_) + (?!+- )B(a+ , a_ + 1 ) 

- «) + (2?i_ - ?i)(?io -«-) + (?!+- ?i-)(«o - »-) (a|+'a_+i) 
«o — n 

One can easily show that n = no in the limit Xi —>■ 00, and « = + for X\ = 0. Finally one 
obtains 

X 

n = nQ + {p+n++p-n- -nQ)— — , (52) 

T + T 

where x = Xj^' and x = ~ 7^ + 77- Thus, n depends monotonously on x. The expression (52) as well 
as the formula for po{n), (^) have been mentioned without a detailed derivation in |^. 
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3.2 The general two-state model 

For a general noise process z{t) the effective rates are 

n = I dz{l - r(z)M-i)-V,-(z)(^o(z), (53) 
where again ri := r2i = — rn, := r\2 = —^n- This yields directly 

__ J^z(l-r(z)M-')-V2(zHo(z) 
/c/z(l-r(z)M-')-V(zHo(z)' 

Whether or not this expression can be calculated explicitely depends on the generator of the noise 
process z.(t). A case where the calculation is possible is a kangaroo process. A kangaroo process is a 



process where all non- vanishing eigenvalues of My are equal. In that case (54) simplifies to 



__ /Jz(l+xr(z))-V2(z)^o(z) 

/Jz(l+Tr(z))-ir(z)(t)o(z) • ^ ' 

The derivation of the last result is not straight forward: One cannot simply replace by — T"^ But it 
turns out, that the result is the same, since additional factors in the denominator and in the numerator 



cancel each other. ( p5[ ) is valid for a dichotomous process as well and yields directly (|52|). 



3.3 Comparison with the Fokker-Planck equation 

As already mentioned, the validity of rate equations for fluctuating potentials is less clear than for fixed 
potentials. It is therefore important to compare results for rate equations with exact solutions of the 
Fokker-Planck equation. They can be obtained easily for one-dimensional, piecewise linear potentials. 



This has been shown in detail for periodic potentials in [12, 13]. In the present case the procedure is 
similar, let us therefore describe it only briefly. For this comparison, we restrict ourselves to the stationary 
distribution /?o(x) = / fifzp(x,z) and to the mean escape times for the particle sitting in a minimum of the 
potential. 

We first divide the domain Q., on which the Fokker-Planck equation is defined, into a set of intervals 
/(, / = 1, . . . ,A^/. We assume that the force /(x,z) = —^^-57^ is constant on each interval //, 

f(x,z)=m ifxG/,-. (56) 

As a second step, we expand the stationary solution of the Fokker-Planck equation using the eigenbasis 
ofM^, 

p(x,z) = A)(x)^o(z) + £ ft W(^/t(z). (57) 

*:>0 

On the interval /,, the Fokker-Planck equation yields 

/;-(z);^;, (x)^o(z) + Ti^l,{x)Uz) + 1 r^l{x)m^^{z) 

k>0 

+T'£p'l'{x),^u{z)-'£hp'k{^)<^k{z) = 0. (58) 

k>0 k>0 

This equation can be integrated once with respect to x. Furthermore, expression of the form fi{z)(^kiz) 
can be expanded, 

fi{zmz)=Y,fi'l,<^k'{z). (59) 

k' 
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Figure 1: The mean occupancy in the left minimum of a piecewise linear, dichotomously fluctuating 
potential as a function of x for various T . The parameters are T = 0.2, 0.4, 0.6, 0.8, 1.0, the average 
force takes the values 10, 3/2, -3/2, 4/3, -4, -10, A/ takes the values 0, 1/4, -1/4, 2/3, -2, 0. The force 
jumps at the x = —4, —2, 0, 3, 4. The maximal value of the mean occupancy decreases with increasing 
T . The dashed lines are the results from the effective rate equation. 

One then obtains 

f&Poi^) + I + TM = 0, (60) 

/t'X) 

/oiVo(x) + ^ f^^p/^ (x) + Tpl{x) = hPkix) , k>0. (61) 

k'>Q 

The coefficients fj^'l, and therefore the final set of equations for pk{x) do not depend on the representation 
of the noise process z{t) as ist should be. If the space of functions ^k{z.) has finite dimension N, ( |60| ) and 
(^ is a set of A'^ differential equation with constant coefficients, which can be solved explicitely. One 
obtains a simple eigenvalue problem for an N xN matrix, which can be solved numerically. A general 
solution of this problem is a linear combination of the N different solutions of the eigenvalue problem. 
The remaining problem is thus to determine the NjN coefficients in these linear combinations. They are 
determined by the normalization of po{x), by the continuity of po{x) and by the continuity of pk{x) and 
p'i^{x) for ^ > 0. This is a problem of solving A'^^A^ linear equations for NjN unknown coefficients, which 
can again be done easily numerically as long as NjN is not too large. In Fig. 1 we show explicit results 
for a one-dimensional dichotomous two-state model. The potential V{x,z) = {V){x) +zAV{x) has two 
minima. z{t) is a dichotomous process which takes the two values ±1. In this case N = 2 and Nj = 6. 
The plot shows the probability to find the particle in the left potential well as a function of x. It can be 
compared to n, which can be calculated from the effective rate equations. A similar comparison for a 
different, dichotomously fluctuating potential was shown in [^. The results from the rate equation yield 
an accurate approximation for sufficiently large x, as expected. The intra-well relaxation time can be 
estimated from the time, the particle needs to reach a minimum of the potential if it starts close to the 
maximum. In our units it is somewhat larger than unity. The approximation becomes better for smaller 
temperature. The behaviour is similar for other dichotomously fluctuating potentials with two minima. 
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X 



Figure 2: The mean occupancy in the left minimum for a fluctuating potential as in Fig. 1. / and A/ 
are the same as in Fig. I, T = 0.4. The different curves are for sums of A^ — 1 dichotomous processes, 
N = 2, . . . ,6. For increasing N the maximum of the curve is shifted to the left. The dashed lines are the 
results from the effective rate equation, the results for A'^ = 4, 5, 6 cannot be distinguished. 



For other stochastic processes z{t) one can draw similar conclusions. Let us consider as an example 
a potential V{x,z) = {V){x) +zAV{x) where now z{t) is a sum of N — I dichotomous processes, where 
each dichotomous process takes the values ibl/^//V. This means that z{t) takes the values {—N + 2n + 

1)/ \//V, n = 0, . . . ,N with the probability p„ = 2^^ ^ ^ ^ . In the limit N ^ oo one obtains an Ornstein- 

Uhlenbeck process [14|. For {V) and AV we take the same values as in Fig. 1. Results for the mean 
occupancy are show in Fig. 2 together with the corresponding results from the rate equation, obtained 
from (54). As for the dichotomous process, the results from the rate equation agree quite well with the 
results from the Fokker-Planck equation for not too small values of x. Similar results can be obtained for 
other noise processes as well. For kangaroo processes, the mean occupancy has been obtained explicitely 
(^. The corresponding results for the mean occupancy are similar to what has been obtained so far. If 
one takes for instance kangaroo processes with the same <^o{z) as for the sums of dichotomous processes 
in Fig. 2, the quantitative results are almost the same, they differ by less than 1%. 

The results in Fig. 2 show that already for small N one obtains essentially the result for the Ornstein- 
Uhlenbeck process. The convergency to the Ornstein-Uhlenbeck process is very fast. For N > 6 the 
curves lie on top of the curve for N = 6. Clearly, the convergency depends on the parameters of the 
system. For smaller temperatures we observe a slower convergency. 

It is possible to consider potentials with more than two minima as well. The corresponding rate 
equations describe systems with more than two states. For simple systems with few minima and e.g. a 
dichotomously fluctuating potential, one can again compare the results from the Fokker-Planck equation 
with the results from the rate equations. The generic behaviour of such systems is again described by the 
effective rates as long as x is not too small. Whether or not ( |29| ) can be used to calculate the effective rates 
analytically depends on the structure of the problem. For complicated noise processes or for potentials 
with many minima, /?eff can only be obtained numerically. 
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4 Mathematical Aspects 



In the main part of this work we investigate the stationary solution of a Fokker-Planck equation (Q) with 
a fluctuating potential, assuming that such a stationary solution exists. The aim of this section is to proof 
the existence of a stationary solution under certain assumptions that will be sufficient for our purpose. 
For the special case of a one-dimensional, dichotomously fluctuating potential, Pechukas and Ankerhold 
[ [lO| ] proved the existence of a stationary solution under some special assumptions on the potential. The 
result of this section can be understood as a generalization of their results. 



4.1 The discrete case 



Let us first discuss the discrete case, which is very simple but shows some interesting aspects. Let us 
assume that a system has A'^ states and that it can move from one state to another. The behaviour is 
described by a rate equation 

^ = JlrijPj. (62) 

We now assume that the rates fluctuate between A^^ different states a = 1 , . . . ,Nr. The fluctuation between 
these states is determined by a matrix M = (moip)a.(3=i that contains the rates for the fluctuations. 
The probability of the fluctuating system to be in the state given by / and a is determined by 



dph 



dt 



m, 



(63) 



P 



The stationary solution of this equation, if it exists, must be the eigenvector of the matrix 



/ Ri+m\il m\2l m\j,I 
m2\I R2 + m22l ni23,I 



R 



msil 



mil Ri + 



mwJ \ 



(64) 



to eigenvalue 0. Here, / is the N xN unit matrix. The matrices /?« are rate-matrices. As a consequence. 



- 52 rjia 



(65) 



and rija > for / ^ j. Furthermore, M is a rate matrix, niaa = — Lp^a'^p.a m^ p > for a 7^ p. 
Therefore, all the off-diagonal matrix elements of R are non-negative. I assume that starting from some 
state of the system, any other state can be reached dynamically. This means that R is irreducible. Then, 
one can introduce a constant c that is larger than the modulus of any diagonal matrix element of R. Let 
/ be the NrN x NrN unit matrix. Then the matrix R + cI is irreducible and has only positive matrix 
elements. As a consequence, the Perron-Frobenius theorem applies. The eigenvalue with the largest 
modulus of ^ + c/ is non-degenerate, real, and positive, and the corresponding left and right eigenvectors 
have only non-negative entries. Now, the vector (1, 1, . . . , 1) is a left eigenvector of R with the eigenvalue 
c. Therefore c is the eigenvalue with the largest modulus of ^ + c/. This shows, that is the eigenvalue of 
R with the largest real part and that its right eigenvalue has only non-negative entries. It is the stationary 



solution of (|63|). One can show that any solution of the rate equation ( p3[ ) tends to the stationary solution 
for t 00. 
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This proof is very simple, it is based on the fact that rates are non-negative and on the Perron- 
Frobenius theorem. Let us now take a look at the original Fokker-Planck equation (|5|). For simplicity, we 
discuss only the case where the potential fluctuates between Nr different states, z takes Nr different val- 
ues. Let us assume that the Fokker-Planck equation is defined on a finite open domain Q. with reflecting 
boundary conditions on the boundary dQ.. If one introduces some kind of coarse-graining, i.e. a partition 
of n into small parts H,, it is then possible to apply the above result and a stationary solution exists. 
This holds for any partition of D. into small parts. Therefore one should expect that a stationary solution 
of the Fokker-Planck equation exists as well. We will show this in the next subsection using a different 
approach. One has to show that the Fokker-Planck operator has a unique right eigenstate with eigenvalue 
and that any solution of the Fokker-Planck equations tends to that solution for t —>^ oo. Since we are 
dealing with a differential operator and not with a matrix, it is not possible to use a Perron-Frobenius type 
argument. 



4.2 The continuous case 

Let us first show that the ratio of any two solutions of the Fokker-Planck equation tends to unity in the 
limit t —>■ oo. As for the discrete case we assume that z takes Nr values Za, OL= I,... ,Nr and we let 

p{x,Za,t) = Pai^j). (66) 

Furthermore, we again assume that the Fokker-Planck equation (^ is defined on a finite open domain Q. 
with reflecting boundary conditions on the boundary 3n. The Fokker-Planck equation can be written in 
the form 

— ^ =V{VVa{x) + TV)pa+Y.'^o^?'P?' = ^^Po'~^ll'^a?,Pf,- (67) 

p p 

Let p\a, P2a be two positive solutions of this equation. We define 

H{t) = £ / d'^xpyaHPla/pia). (68) 



This function is often called Kullback information and has been introduced by Kullback [15, 16]. By 
standard methods one can show that (i) H{t) > and (ii) ^ < 0. As a consequence H tends to zero for 
t oo and therefore pia/pia tends to 1. The proof is similar to the one presented in chapter 6 of [17|. 
First, one shows 

H{t) = £ / d''xp2a f ^In ^ - ^ + l) > 0, (69) 
a-' \P2a P2a P2a J 

which follows from rlnr — r + 1 > for r > 0. Now, one can calculate 
dH f ,d tdpia, t I \ t I \dP2a. 



= Y.Jd''x{-^\n{pia/P2a)-{Pla/P2ay- 



dt 

Pladp2a 
P2a dt 



d'^xl (L(xPia)ln(7?i(x/p2a)+Xl'"ospPpln(pia/P2a 
a \ p 

£ J d'^X fpia(^aln(;?ia/;?2a)+£niap^'pln(7?la/7?2a)- j • (70) 



The first term may be evaluated using 



2 

Llln{pya/P2a) = —V^{pia/ P2a.) " T^i^ {pxa/ P2a)f ■ (71) 
P\a Pta 
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This yields 

^ = -r£ / d'xp,aiyHp,a/p2a)f+ [ d'xl^m^p fp In^^-p ElA . (72) 
dt a-J J ^ P2a P2aJ 

The first term is clearly negative or 0. To investigate the second term, let us introduce 

Aap = §ap+e'Wap, (73) 

where 8 should be sufficiently small (8^^ > maxa Imaal, I will take the limit 8^0 below). Let pia = 
Lp^apPip- One can show that 

/ d'^x)^ pia^n Plain— = / ^^;c> Aappipln \ 

J TV Pia PlaJ J % Pl^Pla 

> /^^.lA„pPipfl-^)=0. (74) 
Expanding the left hand side in powers of 8, one obtains 




and therefore ^ < 0. The limit Pa{x,t o°)of a positive solution of the Fokker-Planck equation is 
positive. If the potential Va{x) is bounded within the domain Q., and if Q. is connected, this Umit yields 
the unique stationary solution of the FoKker-Planck equation, since the Fokker-Planck operator does not 
depend on time. This shows the existence of the stationary solution of the Fokker-Planck equation. 



5 Outlook and conclusions 

One main result of this paper is that the over-damped motion of a Brownian particle in a fluctuating 
potential can be described by kinetic rate equations if (i) the potential fluctuations are slower than the 
relaxation of the particle within a minimum of the potential and if (ii) the positions of the minima of the 
potential do not fluctuate. For temperatures small compared to typical barrier heights of the potential, the 
quantitative agreement of the two different descriptions is very good if one calculates stationary or quasi- 
stationary properties. This shows that the long-time behaviour of a Brownian particle in a fluctuating 
potential has universal properties and does not depend on the details of the potential but only one the 
rates for the transition over the various fluctuating barriers. This is important, since in many realistic 
situations details of the potential are not known. For instance, in the case of a cell surface receptor or 
some other protein in a cell membrane, one knows eventually something about the stable or metastable 
conformations of the protein, but (almost) nothing about the potential that describes the energy of the 
deformations of the protein. Even if the potential is known, it is often simpler to solve a kinetic rate 
equation instead of a Fokker-Planck equation. 

Starting from the kinetic rate equation with fluctuating rates, we derived a formula for effective rates 
which holds for general noise processes. This formula yields directly information on the stationary and 
long-time properties of the system. It can be evaluated for many noise processes. The special case of 
a dichotomous two-state system can be solved exactly and one obtains the distribution functions for the 
occupancy of the two states. 
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As a byproduct, we could show the existence of a stationary solution of the Fokker-Planck equa- 
tion with a fluctuating potential under some general conditions and thereby generalize a result published 



recently by Pechukas and Ankerhold []10|]. 

As a direct consequence of our calculations, one can investigate where rate equations can yield good 
results and where not. This is important since in many applications, mainly in biological situations, 
rate equations are used for a direct comparison with experiments. For instance, dynamical properties of 
membrane proteins are often described and investigated by rate equations. An interesting and well known 
example is the effect of time dependent electric fields on such proteins. It has been discussed it detail by 



Astumian and Robertson [ ]18[ ] using rate equations. In such systems one observes strong amplifications of 
weak signals. Moss [[l^] suggested that the reason may be stochastic resonance. Kruglikov and Dertinger 
[|^] gave a qualitative discussion supporting the occurence of stochastic resonance in such systems using 
a time-dependent potential, but a description using rate equations would have been possible as well. As 



already mentioned, Petracchi et al [ ]11[ ] studied the effect of time dependent electric fields on membrane 
proteins experimentally. As a specific example they used a K+ channel. They measured various stationary 
probabilities and dynamical quantities to describe the statistics of the transitions in these systems. As an 
interesting result they found a phase anticipation: One of the transitions occurs with a negative phase 
shift compared to the stimulus. They compared their experimental results with numerical results for a 
simple two-state Markov model, described by a rate equation. The rate equations do not show the phase 
anticipation. The authors argued that this effect has a biological origin and discuss various possible 
hypotheses to explain it. 

From our discussion it is clear that rate equations do describe the stationary or long-time behaviour 
of a system quite well. Indeed, the stationary or quasi-stationary properties calculated by Petracchi et al 
[ pl| ] agree very well with their experimental findings. The phase anticipation is a dynamical effect that 
occurs on shorter time scales. Therefore one cannot expect that it can be obtained using rate equations. 
Our calculations suggest that this effect depends on the details of the potential and not only on the rates, 
which are determined mainly by the barrier heights. The main problem is that in the case of a specific 
membrane protein the potential that describes the dynamics is not known. Therefore it is not possible to 
model the K+ channel using a time-dependent potential. But it would be of general interest to investigate 
whether or not and under which conditions a Brownian particle in a time-dependent potential shows the 
described phase anticipation. This is clearly beyond the scope of the present paper. 
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